The c-axis charge travehng wave in coupled system of Josephson junctions 
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We demonstrate a manifestation of the charge traveling wave along the c-axis (TW) in current volt- 
age characteristics of coupled Josephson junctions in high-T^ superconductors. The branches related 
to the TW with different wavelengths are found for the stacks with different number of Josephson 
junctions at different values of system's parameters. Transitions between the TW branches and the 
outermost branch are observed. Time dependence of the electric charge in the superconducting lay- 
ers and charge-charge correlation functions for TW and outermost branches show different behavior 
with bias current. We propose an experimental testing of the TW by microwave irradiation. 



Many physical and biological objects may be consid- 
ered as coupled nonlinear oscillators, jljj Particularly, in- 
trinsic Josephson junctions (JJ) in HTSC are studied 
from this point of view often. [2] An attractive feature 
of intrinsic JJ is a generation of coherent radiation in 
terahertz region and it is investigated intensively. 011 A 
coupling between junctions plays an important role. It 
leads to multiple branch structure in the current volt- 
age characteristics (CVC) with additional to single junc- 
tion parameters like the breakpoint (BP) current, the 
BP region width and the transition currents between 
branches. [5- 8] It was shown that the BP in CVC is a 
manifestation of the parametric resonance in this sys- 
tem when the Josephson oscillations create a standing 
longitudinal plasma wave perpendicular to the supercon- 
ducting layers. The parametric resonance features were 
observed experimentally recently. !9'] 

A well known feature of single long Josephson junc- 
tion is the traveling wave along the barrier. 2j It was 
noted 1^, 11 1 that the system of capacitively coupled J J 
has a plane wave which corresponds to the longitudi- 
nal Josephson plasma propagating perpendicular to the 
junctions. The sharp resonance peak in the microwave 
ab sorp tion spectrum obtained in Bi-2212 by Matsuda et 
al.[12| has been identified with this collective Josephson 
plasma mode. Instability of the longitudinal plasma os- 
cillations was discussed in Ref.[llJ and it was proposed 
that the multiple branches are caused by the freezing of 
the collective LPW into the static charge density wave. 

In this letter we demonstrate the additional branch 
in CVC which corresponds to the charge traveling wave 
along c-axis (TW) when the junctions are in rotating 
states (phase differences are not small). Such branches 
were not investigated before. Transitions in hysteresis 
region from the outermost branch (OB) to the traveling 
wave branch (TWB) and between the TWB with dif- 
ferent wavelengths are found. We demonstrate the ap- 
pearance of the TWB dX I — Ic as well which originates 
from different branches of CVC including the zero volt- 
age branch. A different dynamics of correlation functions 
for the TWB and for the OB is shown. We find the ef- 
fect of the microwave radiation on TWB, which allows 



the experimental testing of this branch. 

To simulate the CVC we investigate the phase dy- 
namics in the framework of the capacitively coupled 
Josephson junctions with diffusion current (CCJJ-I-DC) 
model. [13I Il4 | which is determined by the system of 
dynamical equations^ V; = / — sin</j; — Z?^, ^<^/ — 
V; — a(T^+i +Vi-i—2Vi) for the gauge-invariant phase dif- 
ferences ^i{t) between superconducting layers (5- layers). 
Time t is normalized to the inverse plasma frequency uj~ ^ 
(cjp — 2eIc/hC), the voltage - to the value Vq = hujp/2e, 
the current - to the critical current I = Ic- We con- 
sider the periodic boundary conditions in this paper. 
The structure of the simulated CVC in the CCJJ-f DC 
model is equidistant in agreement with the experimen- 
tal results. [15| The details of simulation procedure are 
presented in Refs[8, 16 1 

The manifestation of the TWB is shown in Fig. [1] 
where we present the total branch structure of CVC for 
the stack with 10 coupled JJata — 1, (3 = 0.2. It 
was obtained by multiple sweeping of the bias current 
through the stack. 0, S [i3l This CVC is characterized 
by switching all junctions to the rotating (R-) state at 
/ = Ic, large hysteresis, and multiple branching in the 
hysteresis region by switching the JJ to the oscillating 
(0-) state with the decrease of the bias current. 

The inset (a) enlarges the part of CVC corresponding 
to the transition from the OB to the another branch. As 
we can see, there is one more branch additionally to the 
OB. We prove below that this is TWB. Such TWB was 
manifested in CVC of the stack with 10 J J in Ref.Q (see 
Fig. 5), but it was not mentioned at that time. 

To see the origin of TWB we investigate a transition 
region in CVC near Ic. In during of the sweeping bias 
current procedure we observe the jumps to the OB and 
TWB both from the zero voltage state directly and from 
the other branches of CVC. The inset (b) enlarges the 
part of CVC corresponding to this transition region. 

The TWB have been observed in the stacks with differ- 
ent number of junctions. Particularly, the TWB in CVC 
were found in the stacks with 14 J J and 28 J J (Fig. [2]), 9 
J J (Fig. [3ji) and 7 J J (Fig. SJ. In case of even number of 
junctions in the stack, the jump from the OB to the TWB 
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FIG. 1: (Color online) Manifestation of the TWB in CVC of 
the stacks with 10 JJ. The insets enlarges the parts of CVC: 
(a) near transition from the OB to another branch; (b) near 
transition at I — Ic. where the TWB is created. 



happens first to some unstable state. Such transition for 
stacks with 14 J J at a — 1, (3 — 0.2 is demonstrated in the 
inset (a) of Fig. [21 where a transition from the parametric 
resonance region of OB is observed. For the stacks with 
odd number of junctions we observe usually the jumps 
to the TWB from the chaotic part of breakpoint region 
in OB. The increase of bias current after jump instead of 
its decrease allows us restore the total TWB (see Fig. [3]) . 




FIG. 2: (Color online) Demonstration of the TWB corre- 
sponding to the traveling waves with different wavelength, 
(a) The manifestation of the TWB for the stack with 14 JJ. 
Inset demonstrates the transition from the OB to the first 
TWB with A = UD. (b) The same for the stack with 28 JJ. 

(c) The distribution of the charge in S-layers along the stack 
with 28 JJ for TWB with A = Ao at bias current I = 0.577. 

(d) The same for A = Ao/2. (e) The same for A = Ao/4. Fig- 
ures (f)-(h) show the corresponding Lissajous charge-charge 
diagrams. 



We note that the longest wavelength of a TW in the 
system with N J J is Aq — ND, where D is a period of 
lattice ( the sum of thicknesses of insulating and super- 
conducting layers). So, the waves with A = ND/n might 
appear in the stack, where n is devisor of number N. We 
stress that the TWB in Fig.[T]i(;as naturally obtained just 
by sweeping of bias current. Another possibility is to use 
the variation of the initial conditions for phase differ- 
ences. In Fig. [2] we show such realization of the TWB 
corresponding to TW with different wavelengths. Par- 
ticularly, in Fig. [2^ we demonstrate the manifestation of 
the TWB in CVC of the stacks with 14 JJ related to the 
waves with Aq = 141? and A = Ao/2. Inset demonstrates 
a transition of the OB to the first TWB with A = Aq. The 
analysis of the time dependence of the charge in S-layers 
and its distribution along the stack around the second 
transition in Fig. [2^ shows that the TW with A = Aq 
transforms to TW with A = Ao/2. 

In Fig. EJb) we demonstrate the TWB for the stack 
with N=28, corresponding to the states with wave- 
lengthes Ao = 281?, Ao/2 and Ao/4 . The different initial 
conditions for the phases and the voltages in the JJ were 
used to obtain these branches. We find the wavelength of 
the TW by the distribution of the charge in the S-layers 
along the stack at fixed time moment. Such distribution 
at bias current / = 0.577 is shown in Fig. [2rc),(d) and 



(e). These figures demonstrate also that the amplitude 
of charge is greater for the TW with smaller wavelengths. 

To test if these distributions follow precisely to the 
TW, we present in Fig. [2][f),(g) and (h) the correspond- 
ing Lissajous charge-charge diagrams. In fact, two dif- 
ferent plots are shown in each figure: squares show the 
{Qi, Qi-i) for different values of I at fixed time moment; 
the curves show the trajectory for fixed value of I in dur- 
ing a whole time domain (i.e., the time domain at con- 
stant value of current). We got the same curve for differ- 
ent I, which means that the charge in each layer has the 
same frequency and amplitude, but the charge in layers 
differers by a phase shift only. So, this wave is not a 
standing wave, but a traveling one. 

An example of the charge dynamics corresponding to 
the OB and TWB for the stack with 9 JJ at a = 3 and 
/3 = 0.2 and its simulation procedure is demonstrated in 
Fig. |3l First we follow the OB (arrow I), decreasing the 
current from the value 1=1.2. At 1=0.78405 a transition 
from the chaotic part of OB to TWB is happened. We 
continue the current decreasing process along TWB till 
the 1=0.78. At this point we increase current (arrow 
2), follow the TWB till / = 5. Then we return back 
(arrow 3) and follow along the TWB till the jump to 
the another branch (arrow 4). The insets (a,b) show a 
time dependence of charge in first (solid line) and second 
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FIG. 3: (Color online) Manifestation of the TWB in CVC of 
the stack with 9 JJ at a = 3, /3 = 0.2 and periodic BC. The 
insets (a) and (b) show the charge oscillations in two S-layers 
at 1=0.792 in OB and TWB, relatively. The inset (c) presents 
the results of FFT analysis of charge oscillations for TWB at 
1=0.792. 



(dashed line) S-layers of the stack at / = 0.792 for the 
OB and TWB, correspondingly. We see that in the case 
of TWB the charge oscillations in S-layers are just shifted 
in phase. In the supplement we present the animation of 
the charge-time dependence demonstrating the traveling 
wave along the stack. The inset (c) presents the result of 
FFT analysis of charge-time dependence corresponding 
to the TWB. The time domain was taken at 1=0.792, 
where the value of voltage is equal to V=32.956. As 
we see the value V/N coincide with the results of FFT 
analysis ^ = 3.912. 

One more manifestation of the TWB is shown in Fig.|4l 
where we demonstrate a transition to the TWB for the 
stack with 7 JJ at a = l,/3 = 0.2. The charge-time 
dependence demonstrates a changing of character of the 
charge oscillations in S-layer near a transition point at 
1=0.5509 (Fig. |41d). We can see the transformation of 
the chaotic behavior in the OB to the regular one in TW 
state. In Fig. lit we show the Lissajous charge-charge di- 
agram in chaotic (thin line) and TW (thick line)regions. 
The diagram demonstrate the variation in time of the 
charges in two neighbor layers. The charge in the first 
S-layer Qi is plotted along the x-axis and charge in the 
second one Q2 along the y-axis. As we can see, the Lis- 
sajous charge-charge diagram presents an open trajectory 
for the chaotic region and closed one for the TW region. 

We find the effect of microwave radiation on TWB, 
which was taken into account by the including an addi- 
tional current Ir = A* sinw^i through the stack. In- 
creasing power of radiation ( amplitude A ) leads to the 
destroying of the TW and transition from the TWB to 
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FIG. 4: (Color online) (a) Manifestation of transition OB — >■ 
TWB in the CVC of the stack with 7 J J. (b) The charge-time 
dependence for S-layer at transition point, (c) Charge-charge 
diagram for the neighbor S-layers in chaotic region of the OB 
and in TWB. (d) Effect of microwave radiation on the TWB. 
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FIG. 5: (Color online) The dependence of correlation function 
"charge-charge" for neighbor layers three and four on the 
bias current. The procedure of simulation is the same as for 
Fig. [3] Inset (a) stress the growth of C for the TWB in the 
bigger interval of I. Inset (b) enlarged the C'^{I) dependence 
near the transition from the OB to the TWB. The thick curve 
show the CVC in this transition region. 



the another branch in CVC. In Fig. |3Ji we demonstrate 
the effect of radiation, which was switch on at J = 0.55 
with amplitude A ^ 1 and Josephson frequency loj. As 
we can see, the TWB disappears under this radiation. 
This result open a way for the experimental testing the 
observed branches in CVC. 

To compare dynamical behavior of the states cor- 
responding to the OB and TWB in long bias cur- 
rent interval, we have made a correlation func- 
tions analysis [1^. We have investigated a correla- 
tion of charge in neighboring layers by charge-charge 
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FIG. 6: (Color online) The charge-time dependence in transi- 
tion region near I = Ic (see the inset (b) in Fig. [1} . The inset 
(a) enlarges the Qiit) dependence in three first S-layers in 
zero- voltage state. The inset (b) demonstrates the Qiit) de- 
pendence for the same S-layers in the transition region shown 
by oval. The inset (c) shows the Qi(t) dependence for the 
same layers in the state with TW. 
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=< Qi{t)Qi+i{t) >= 
Qi(t)Qi+i{t)dt. Fig. E] shows 
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the dependence of correlation function "charge-charge" 
C"^ for neighbor layers three and four C| 3 on bias cur- 
rent for the OB and the TWB. We stress that all curves 
were obtained in the one numerical experiment. Inset 
(a) demonstrate the increase of correlations (C^) in the 
TW state with decreasing of the bias current in the 
large interval of /. As it has been seen from the in- 
set (b) which enlarges the region near BP, the corre- 
lation functions C"^ for the OB demonstrate the char- 
acteristic BPR structure with the features obtained in 
the Ref. 16]. We see also that the features of correla- 
tion functions coincide with the features of CVC. The 
Cii_i for TWB are coincide for all /. The correlation 
functions current-current (not presented here) show the 
same behavior. The autocorrelation functions for TW 
state, which are defined by Cf =< Qi{t)Qi{t — ti) >= 

limTm-n^oo Tm-T» " ti)d't: have a con- 

stant amplitude, i.e., the system in the state with TW is 
complectly self repeating. 

Finally, let us describe shortly the creation of the TWB 
at / = /c. In Fig. [5] we present the charge-time depen- 
dence in transition region around I = Ic (see the inset (b) 
in Fig. [T]) for three S-layers in the stack with 10 JJ. We 
enlarged the Qi (t) dependence in this region and show it 
in three insets separately: the inset (a) shows it in zero- 
voltage state, inset (b) - in some part of the transition 
region shown by oval, inset (c) - in the state with TW. 
The simulation was done by increase a bias current in the 
interval (0.993, 1.003) around I = Ic with a step in cur- 



rent 51 = 10 We observe here a complex dynamical 
behavior in a very short time interval inside of the time 
domain. Fluctuations of charge in S-layers in zero volt- 
age state is replaced by specific changes of the charge and 
the oscillations with nonzero time averaged value (charge 
value in S-layer is high enough, up to \Q\ = 2.5), then 
we see a transition to the TW state. The amplitude of 
charge oscillations in TW state is large, but its average 
in time is zero. A detailed description of the switching 
dynamics will be presented somewhere else. 

In summary, we showed that CVC of coupled Joseph- 
son junctions has the branches related to the charge trav- 
eling waves along the c-axis. Transitions between such 
branches demonstrating the changing of wavelength of 
the traveling wave with a decrease of bias current were 
found. We fixed also the transitions from the outermost 
branch to the traveling wave branch. Detailed analy- 
sis of the time dependence of the charge in supercon- 
ducting layer, its FFT analysis and the investigations 
of correlation functions showed different features of such 
transitions, particularly, different behavior for the out- 
ermost and traveling wave branches. These results shed 
light on the resonance features of the coupled Josephson 
junctions which are investigated intensively today. The 
detailed analysis of the branch structure in the experi- 
mental CVC of IJJ has not been done yet. We proposed 
a method to distinguish the branch related to the charge 
traveling wave along c-axis from the other branches in 
CVC by microwave irradiation. 

We thank R. Kleiner, A. Yurgens, V. Krasnov, M. 
Suzuki, F. Mahfouzi and M. R. Kolahchi for helpful dis- 
cussions. 
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